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Abstract 

A competitive (2+l)-dimensional model of deposit formation, based on the com- 
bination of random sequential absorption deposition (RSAD), ballistic deposition 
(BD) and random deposition (RD) models, is proposed. This model was named as 
RSADi_ s (RDjBDi_f ) s . It allows to consider different cases of interparticle interac- 
tions from complete repulsion between near-neighbors in the RSAD model (s = 0) 
to sticking interactions in the BD model (s = l,f = 0) or absence of interactions 
in the RD model (s = 1, / = 0). The ideal checkerboard ordered structure was 
observed for the pure RSAD model (s = 0) in the limit of h — > oo. Defects in the 
ordered structure were observed at small h. The density of deposit p versus sys- 
tem size L dependencies were investigated and the scaling parameters and values of 
Poo = p(L = co) were determined. Dependencies of p versus parameters of the com- 
petitive model s and / were studied. We observed the anomalous behaviour of the 
deposit density with change of the inter-particle repulsion, which goes through 
minimum on change of the parameter s. For pure RSAD model, the concentration of 
defects decreases with h increase in accordance with the critical law p oc h~ XRSAD , 
where xrsad ~ 0.119 ± 0.04. 
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1 Introduction 



Recently the study of the nonequilibrium processes of deposit formation be- 
came one of the most active research areas in physics and chemistry because of 
both fundamental and practical interest [1]. The essential progress in this field 
was reached owing to application of computer simulations methods for study- 
ing the problems of adsorption, formation of thin films and coating. There 
exists a limited number of computer models, allowing to carry out effective 
calculations and to establish the main scaling laws for formation of deposits 

[2]- 

Among the most popular models for deposit formation simulation are the 
models of random sequential adsorption (RSA), ballistic deposition (BD) and 
random deposition (RD) [2,3]. In these models, the particles get fixed after 
deposition and don't move, so these models describe the growth processes 
far from equilibrium. It is also assumed that particles are rigid and can not 
overlap. In RSA model, particles get fixed at some distance one from another 
[3,4]. In RSA model with the nearest neighbor (NN) exclusion, configurations 
with the particles having nearest neighbors are eliminated. In fact, it means 
existence of some repulsion between particles. In RSA model, all sites in the 
lattice can be filled with equal probability. However, this requirement is not 
fulfilled on deposit formation. In this case, the previously deposited particles 
can screen free sites, located below. So, we named the variant of RSA model 
for deposit formation as RSA deposition model (RSAD). In BD model, the 
particles stick at the point of their first contact. It means existence of the 
short-range attraction. In RD model, the particles deposit without sticking 
and it means existence of the short-range repulsion. Recently, a number of 
mixed or competitive models were proposed. They are based on consideration 
of deposition for different kinds of particles [5,6,7] 

In this work, the competitive (2+l)-dimensional model with three kinds of 
particles, depositing according to the rules of random sequential absorption 
deposition (RSAD), ballistic deposition (BD) and random deposition (RD), 
is investigated. This model allows to explore a wide class of systems with 
different interparticle interactions, which vary from a complete repulsion as 
for RSAD model to short-range attraction as for BD model. 

The paper is organized as follow. The model is described in section 2. In 
section 3, the scaling behaviour of the deposit density and concentration of 
the structure defects are discussed for different values of the competitive model 
parameters. Concluding remarks are presented in section 4 
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2 Model 



The competitive (2+l)-dimensional model, which combines models of the ran- 
dom sequential adsorption deposition (RSAD), ballistic deposition (BD) and 
random deposition (RD), named by as RSAD i_ s (RD/BDi_/) s , where s and 
/ parameters are fractions of particles with different kinds of interparticle 
short-range interaction potentials. A particle of RSAD or BD or RD kind ran- 
domly falls straight down onto a growing surface, one at a time, and deposits 
at a site of the cubic lattice. 

Parameter s characterizes the short-ranged next-near-neighbor repulsion of 
particles in the deposit, and parameter / characterizes adhesion of particles. 
In the extreme case of s = all the particles are of RSAD kind. This case 
corresponds to the strong repulsion between particles in NN sites and they 
can deposit only in the next NN sites. When s — 1, p — 0, all the particle are 
of BD kind and newcomer sticks to the deposits at a point of its first contact. 
When s — 1, p — 1, all the particle are of RD kind. This case corresponds 
to existence of the short-ranged repulsion between particles and formation of 
completely compact deposits. 

The density of a deposit may be defined as 



where N is the number of deposited particles, and V is the volume of deposit. 

The value of p is calculated at a moment when saturation was reached in the 
system with the volume V = L x L x L. At this moment, filling of the system 
with new particles was terminated. 

For RSAD model, the structure of deposit becomes more regular at large 
deposit height h, and in the extreme case of h — > oo, the ideal checkerboard- 
ordered structure get formed. In this ideal structure, the empty lattice sites 
alternate with occupied ones and each empty (or occupied) site has exactly 
z = 6 occupied (or empty) NN sites. But at small values of h there exist 
some defects in the regular structure (Fig. 1), when NN exclusion number z 
is different from 6. 

The concentration of defects p(h) may be defined as 



where z(h) is the NN exclusion number averaged over all the occupied and 
empty sites in the horizontal layer with the constant h. 



p = N/V 



(1) 



p{h) = l-z(h)/6, 



(2) 
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Fig. 1. Schematic presentation of the deposit for pure RSAD model (s = 0). Black 
cubes correspond to the sites filled with particles. In the ideal checkerboard-ordered 
structure each particle is surrounded by 6 empty sites. The examples of defects are 
also shown. 

The periodic boundary conditions are applied along the horizontal directions 
x and y. The size of the base L was varied within the interval of L = 10 — 400 
and deposit height was h < 200L. The values of p and p(h) were averaged 
over 100 — 1000 of different configurations for each fixed set of parameters s 
and /. 



3 Results and discussion 



At the beginning we studied the scaling behavior of the deposit density p in 
cubes with varying sizes L x L x L. A clear scaling law of the following type 
was observed: 

p=p 00 + aL- 1 '\ (3) 
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Fig. 2. Density of deposit p and scaling exponent v versus s in equation 3 at different 
values of /. Insert shows similar plots of scaling amplitude a p versus s. 

Here, poo = p(L — > oo), v is the scaling exponent and a is the amplitude. 

Figure 2 shows Poo, v and a versus s dependencies at different values of /. As 
we have noted previously, for RSAD model (s = 0) in the limits of L — > oo and 
h — > oo the ideal checkerboard ordered structure get formed, that corresponds 
to poo = 0.5. With increasing of s when / is a constant the density of a deposit 
Poo always goes through the minimum at s = s min . The value of s m j n increases 
with decreasing of /. The point of s = 1 and / = corresponds to a pure BD 
model for which poo = 0.3000 (this value in accordance with data [8]). The 
point of s = 1 and / = 1 corresponds to the pure RD model, when poo = 1 
and compact deposit forms. 

At first glance, the Poo versus s behavior is somewhat anomalous. Really, the 
effective repulsion between particles decreases with increase of s (at constant 
/). So, we should expect compacting of deposit and increase of Poo- Another 
interesting feature is that the scaling exponent v decreases abruptly with s 
decrease in the interval of s < 0.1. In fact, we observe two different scaling 
regimes with large and small values of scaling exponent v. For the pure RSAD 
model (s = 0), scaling exponent is vrsad = 8.4±0.3. In a regime far from pure 
RSAD model, the scaling exponent falls within the interval of v — 1.0 — 1.8. 
The amplitude a continuously grows with s increase and it approaches zero at 
s < 0.1, which results in large errors of v estimation (see insert in Fig. 2). 
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Fig. 3. A typical map of the heights of defect annihilation in the plane of x — y. 
For clearer presentation, a small system of the size 30 x 30 x 1400 was taken. 

The nature of unusual behavior of the deposit density p^ with variation of s 
may be understood more clearly from the analysis of defects evolution during 
the deposit growth. The initial abrupt decreasing of the deposit density with 
s increase (at s < 0.1) reflects generation of defects in the ideal checkerboard- 
ordered structure. Naturally this results in loosening of the structure of de- 
posit. The point of density minimum at s = s m i n corresponds to a kind of 
equilibrium between the processes of birth and vanishing of defects. With fur- 
ther s increase, effectiveness of regeneration of the ideal checkerboard-ordered 
structure decreases. Finally, it results in increase of p^ with s increase. 

The described behavior of p^ can be compared with behavior of the concentra- 
tion of defects p(z). Figure 3 shows the typical map of the defect annihilation 
heights in the plane of x — y for the pure RSAD model. Here, the lighter color 
corresponds to the higher deposit coordinate z, at which defects disappear in 
the ideal RSAD structure. 

At initial moment of deposit formation, there arise a lot of defects in the ideal 
RSAD structure. This is a result of the absence of correlations between depo- 
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Fig. 4. Concentration of RSAD defects p versus the deposit height h at / = 1 and 
different s and L. The dashed line corresponds to the power equation (4) with the 
slope xrsad ~ 0.119±0.04. The value of corresponds to the limit h — ► oo. Insert 
shows the limiting values of concentration of the RSAD defects (at h — > oo) and 
density of a deposit (at /i = L — > oo) versus s at different values of /. 



sition in the columns at small height of the deposit. These defects can spread 
in the direction of deposit formation along z axis and their concentration in 
each layer p(z) increases or decreases owing to the processes of defect birth 
(at s ^ 0) or annihilation, respectively. 

Figure 4 shows the examples of defects concentration p versus deposit height 
h dependences for the case when / = 1. The similar dependences were also 
observed for other values of /. 

For pure RSAD model (s = 0) the concentration of defects p continuously 
decreases with increase of the deposit height h. In the limit of h — > oo the 
defectless checkerboard-ordered structure get formed. In the latter case, the 
obvious scaling behavior of p(h) dependencies is observed for systems with 
different size of base L. We believe that a simple power law can be applied for 
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description of p(h) behavior in the limit of L — > oo 

p oc h~ XRSAD (4) 

where Xrsad is the scaling exponent of defects annihilation for the pure RSAD 
model. 

This behaviour is rather similar to the critical annihilation of inter-domain 
boundaries in the model of competitive growth by Saito and Muller-Krumbhaar 
(SMK) [9] ( see also, [10,11]). 

Prediction of the scaling exponent in critical annihilation is not trivial. For 
normal Brownian motion of annihilating defects the scaling exponent should 
be xb = 1/2 [12]. In 1 + 1 dimensional SMK-model another scaling exponent 
Xsmk ~ 2/3 [9] was observed. This result was explained by existence of dif- 
ferent mechanisms of the critical annihilation of defects in SMK and random 
walk motion models. It was conjectured that the lateral displacements of inter- 
domain boundaries in the SMK model are controlled by the same processes 
as those causing roughening the outer interface width and Xsmk ~ where 
a is the roughening exponent. On the other hand, the value of a can reflect 
the growth mechanism, geometrical confinement of growing interface, space 
dimension, etc. [13]. For example, for usual models of KPZ class of universality 
[14] the roughening exponent is a — 1/2 in 1 + 1 dimension and is ~ 0.2 — 0.4 
in 2 + 1 dimension [2] . 

For pure RSAD model (s = 0) and an isotropic system (L x L x L), we can 
estimate the value of the scaling exponent in Eq. 4 as Xrsad = 1/vrsad ~ 
0.119 ± 0.04, where vrsad — 8-4 ± 0.3 is the scaling exponent for the density 
of deposit in Eq. 3. Indeed, we believe that the scaling changes in the deposit 
density p and in the concentration of defects p are controlled by the same 
process, and p(L) oc p(L). The dashed line in Fig. 4 corresponds to the slope, 
estimated from scaling of the deposit density. But for anisotropic systems 
L x L x h (h > L) , the concentration of defects decrease more quickly than 
it is predicted by the power law in Eq. 4, and the most obvious deviations are 
observed with increase of the system anisotropy degree at high h (see Fig. 4). 

In a general case of the competitive model RSAD ]_ s (RD^BD^j),,, the pat- 
terns of defect distributions in a deposit can be more complicated. Here, defects 
can birth or can vanish and competition between these processes results in for- 
mation of a certain defected structure with the finite concentration of defects. 
For mixed model at s ^ 0, the critical behavior of (4) type disappears. In 
the limit of indefinitely large h (h — > oo), the concentration of defects remain 
finite and has a nonzero value p^. With s increase, which corresponds to the 
weakening of repulsion between the particles, increase of p^ is observed (See 
insert in Fig. 4), and in the limit of pure RD model, (s — 1, / = 1), = 1, 
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because deposit is compact in RD model and even local checkerboard-ordered 
structure is absent. 



4 Conclusions 

It was demonstrated that implementation of the RSA rules for 3d- growth 
of deposit on a cubic lattice with the nearest-neighbor exclusions results in 
formation of a checkerboard-ordered structure. The initial disorder of this 
structure critically disappears when the deposit grows and, finally, an ideal 
checkerboard-ordered structure gets formed. In a competitive RSAD 1 _ s (RDjBD 1 _j) 
model, the deposit density goes through the minimum with variation of s pa- 
rameter at constant /. Such behavior reflects the influence of competitive 
processes related with the birth and annihilation of defects. 
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